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A typical commercial aircraft trajectory consists of a senes of flight segments. An aircraft switches from one 
segment to another when certain specified variables reach their desired values. Trajectory synthesis for air traf- 
fic control automation must be consistent with practical pilot procedures. We examine capture conditions for 
merging trajectory segments to model commercial aircraft descent in trajectory synthesis. These conditions trans- 
ate into bounds on measurements of atmospheric wind, pressure, and temperature. They also define ranges of 
thrust and drag feasible for a descent trajectory. Capture conditions are derived for the Center-TRACON Au- 
tomation System developed at NASA Ames Research Center for automated air traffic control. Various uses of 
captiire conditions are discussed. A Boeing 727-200 aircraft is used to provide numerical examples of capture 


Nomenclature 

A, B = defined constants 

Bi(h) = upper bound on wind gradient 

Biih) = lower bound on wind gradient 

D = aircraft drag 

g = acceleration because of gravity 

h = altitude 

L = aircraft lift 

M = Mach number 

m = aircraft mass 

P = air pressure 

p = air pressure normalized by sea-level pressure 

B - gas constant 

r - wind magnitude bound 

s = path distance 

T = engine thrust 

V cas = calibrated airspeed 

Vg = ground speed 

V, = true airspeed 

V w = wind speed 

x = state vector in a segment 

y = generic capture variable 

8 = lower capture bound 

€ = capture accuracy 

Ya = air-relative flight-path angle 

Yi = inertial flight-path angle 

k = thrust setting 

<p = aircraft bank angle 

^ = inertial velocity heading 

B = atmospheric temperature 

® = temperature normalized by sea-level value 

0rw = relative wind angle 

( )o = constant values 

( ).v = sea-level values 

( )* = desired values 
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Introduction 

G ROWTH in nationwide air traffic has put severe burdens on the 
current air traffic control systems and created heavy workloads 
for air traffic controllers. Further increase in the throughput of the 
existing air traffic control infrastructure will be greatly enhanced by 
computer automation. 

Researchers at NASA Ames Research Center have developed a 
Center-TRACON Automation System (CTAS) for automated air 
traffic control. 1 CTAS generates landing times and descent advi- 
sories for arrival aircraft over extended terminal areas consisting of 
a Center airspace and a TRACON airspace. Center stands for Air 
Route Traffic Control Center and covers the en route airspace area 
about 200 n miles around an airport. TRACON stands for Terminal 
Radar Approach Control and covers the terminal airspace area about 
40 n miles around an airport. Computed advisories allow aircraft to 
descend efficiently and to meet specified arrival times while avoid- 
ing conflict. CTAS has been shown to be able to decrease controllers’ 
workload and to increase aircraft landing rates. 

CTAS comprises several generic modular tools. The Traffic 
Management Advisor schedules landing sequences and times for 
arrival aircraft. The Descent Advisor provides advisories for de- 
scending aircraft in the Center airspace. The Final Approach Spac- 
ing Tool provides advisories in the TRACON airspace for merging 
traffic onto the final approach. All of these tools rely on the aircraft 
trajectories calculated in the Trajectory Synthesizer module. 

An air traffic control automation system has two essential com- 
ponents: a trajectory synthesizer and a scheduler. The trajectory 
synthesizer calculates aircraft trajectories that are efficient, conflict- 
free, and consistent with air traffic control and pilot procedures. It 
also can vary descent speeds to meet scheduled arrival times. The 
scheduler generates a sequence of landing orders and times based on 
aircraft trajectories predicted by the trajectory synthesizer. There- 
fore, the success of an air traffic control automation system hinges 
on the reliability of trajectory synthesis algorithms. 

Trajectory synthesis for air traffic control automation should be 
consistent with practical pilot procedures. In a realistic commercial 
aircraft descent, a complete trajectory is divided into a series of flight 
segments. Mathematically, each flight segment can be described by 
two constant control variables selected from among engine thrust 
setting, Mach number or calibrated airspeed, and altitude rate or 
inertial flight-path angle. 2 - 3 These segments can be specified by 
controllers via voice communication and can be flown by pilots with 
standard instrumentation. Furthermore, airline specifications often 
combine a number of segments in a specified order to form descent 
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profiles. Three types of descent profiles are currently modeled in 
CTAS: fast, nominal, and slow. Each type of profile covers a range 
of arrival times. 

For trajectory synthesis in air traffic control automation, point- 
mass aircraft equations are adequate. Inputs to trajectory synthesis 
include a series of waypoint locations that define the horizontal 
path, initial aircraft states from radar measurements, aircraft thrust 
and drag models, measurements of atmospheric conditions, and 
constraints because of aircraft performance and air traffic con- 
trol. With given descent and/or cruise speeds, aircraft equations 
of motion are numerically integrated segment by segment to form 
complete trajectories. Then, descent and/or cruise speeds are it- 
erated to achieve specified final times. Therefore, outputs of the 
trajectory synthesis are aircraft trajectories that meet specified final 
times. 

For two flight segments to connect properly, certain variables, 
called capture variables, need to reach their desired values. Specif- 
ically, a capture variable is monitored during the integration of a 
segment. Once the capture variable reaches the desired value within 
a certain tolerance, the capture is declared successful and the inte- 
gration of a new segment is started. Each segment has a specified 
capture variable, depending on its location in a descent profile. 

Clearly, the convergence of trajectory synthesis algorithms de- 
pends on three related elements: 1) numerical integration of eveiy 
flight segment, 2) proper connections of flight segments, and 3) it- 
eration over speeds to meet a specified time of arrival. Stability of 
numerical integration depends on the aircraft equations, the integra- 
tion scheme, and the integration step size. In CTAS, a second-order 
Runge-Kutta integration scheme is used and integration step sizes 
are varied according to the nature of each flight segment. These 
choices achieve a proper balance between accuracy and computa- 
tional speed. The search over speeds to achieve a desired time of 
arrival uses function values only and is based on bracketing the 
solution with an ever-decreasing interval. It is limited to four iter- 
ations. If a desired speed cannot be found within four iterations, a 
different horizontal path can be used and the trajectory synthesizer 
module can be called upon again. In contrast, the trajectory syn- 
thesis algorithm cannot obtain a complete aircraft trajectory if any 
two segments fail to connect in capture. Therefore, successful seg- 
ment connections represent a crucial requirement to convergence of 
trajectory synthesis. 

The connection of two segments in a numerical integration is 
guaranteed if the capture variable in the first segment changes mono- 
tonically. In this paper, conditions required for capture variables to 
behave monotonically are presented for CTAS. Capture require- 
ments are unique to the Descent Advisor trajectory synthesis. These 
conditions translate into bounds on predicted weather conditions, 
including, wind components, pressure, and temperature, as well as 
feasible ranges of thnist minus drag. 

The remainder of the paper states the capture requirement math- 
ematically, presents details of capture conditions for different seg- 
ments and profiles, and uses a Boeing 727 aircraft example to illus- 
trate these capture conditions. A summary of key results is given at 
the end. 


Problem Statement 

Trajectory synthesis for the Descent Advisor calculates aircraft 
trajectories from the current aircraft position to a specified boundary 
point (Metering Fix) between the Center airspace and the TRACON 
airspace, for one aircraft at a time. A computed trajectory reaches 
the crossing point at specified altitude, speed, and time of arrival. 4 
Over the Center airspace, hundreds of these four-dimensional air- 
craft trajectories need to be generated every 12 s, which is the cycle 
of radar sweeps. As a result, trajectory -synthesis algorithms need to 
be computationally efficient. In CTAS, the Center Trajectory Syn- 
thesizer integrates a set of point-mass aircraft equations of motion, 
simplified to speed up computations 5 : 


. T(h, M ; k) - D{h, V„ L) ^ d 

V t = — gy a + -r (V w cos^rw) (1) 

m dr 


tl = YaVt=YiV g 


( 2 ) 



Table 1 

Flight segments in a descent 


No. 

Segment 

type 

Integration 

variables 

Constant 

controls 

Capture 

variables 

1 

Cruise 

s 

h = 0, (M or V cas ) 

s'! 

2 

Descent accel 

s , V„ h 

at, {h or Yi) 

V^.or/t 1 

3 

Constant Mach 

h 

M, (at, k or yi) 

Vis or ft* 

4 

Constant Vcas 

s , h 

Vcas, 0^ Yi ) 

or 

vis, or.? 1 

5 

Level deceleration 

s,V, 

K,'h = 0 


Top of Descent 



1 

Level Decel 

Fig. 1 Fast-type descent trajectory for high initial altitudes. 


i = \f — ( Vy, sin — V w cos dp# Vg (3) 


mV s 


and 


L = mg/ cos <f> (5) 

In the National Oceanic and Atmospheric Administration weather 
models currently used in CTAS, wind components (V^, 0 W ), pres- 
sure p, temperature 0, and actual altitude are measured at a grid of 
60 x 60 km horizontal stations. At each station, they are given at a 
series of pressure altitudes. 

In the above equations, (V M h) describes the vertical path and 
(r, 40 describes the horizontal path. The ground speed V* couples 
the horizontal and vertical paths. The horizontal route is synthesized 
first on the basis of an approximate vertical profile. In particular, the 
heading angle 4> is approximated by a combination of straight lines 
and constant-radius turns at waypoints. 6 

The vertical trajectory is divided into a series of flight segments 2 3 
to be consistent with current pilot procedure. Each segment is de- 
fined by setting two control variables constant. These variables come 
from the following three groups: engine thrust setting, speed (Mach 
or calibrated airspeed), and vertical rate (altitude rate or inertial 
flight-path angle). Table 1 lists five flight segments in a realistic 
descent, and for each segment, variables that need to be integrated, 
capture variables that are discussed below, and control variables 
that are set constant. Arrows next to the capture variables indicate 
whether they increase or decrease in that segment with respect to 
real time. 

Furthermore, a series of segments are combined in a predeter- 
mined order to form descent profiles. There are three types of pro- 
files used in CTAS: fast, nominal, and slow. Figure 1 summarizes 
a fast-descent profile for high initial altitudes. In comparison, a 
fast-descent profile for low initial altitudes does not contain a con- 
stant Mach descent segment. A nominal descent profile does not 
contain a descending-acceleration segment. A slow descent pro- 
file can be obtained from a nominal descent profile by including a 
level-deceleration segment after the cruise segment and removing 
the Mach descent segment. Each profile defines a range of possi- 
ble descent-calibrated airspeed variations. The descent-calibrated 
airspeed is varied within these ranges to achieve a specified arrival 
time. 

A trajectory of fast, nominal, or slow type is numerically inte- 
grated segment by segment. The integration is either forward or 
backward in time, and the order in which the segments are inte- 
grated depends on what information is known. In the fast profile 
shown in Fig. 1, both the initial aircraft states and the desired final 
states at the crossing fix are known. For a certain descent speed, the 
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location of the top of descent needs to be determined so that the 
trajectory meets initial and final conditions. As a result, the fast tra- 
jectory is integrated backward from the level-deceleration segment, 
as indicated by the arrows in Fig. I. Then, the descent-calibrated 
airspeed and/or descent Mach are varied to meet a specified final 
time. Currently in CTAS, the search over descent speed is based on 
bracketing the descent speed with an ever-decreasing interval, and 
is limited to four iterations. 

Vertical trajectory synthesis involves the calculation of time his- 
tories for three variables: true airspeed, altitude, and path distance. 
The specifications of a flight segment often remove the need to inte- 
grate one of the three variables directly. For example, true airspeed 
can be solved algebraically in a constant Mach descent or calibrated 
airspeed descent segment. Mach number and calibrated airspeed, 
V cas , are related to true airspeed as follows 7 : 

V, = MjlARQ(h) (6) 


V L = 7K0,((p(A)[(O.2A/ 2 + l) 35 - 1] + l}v - l) (7) 



In integrating a segment, a certain variable is monitored to de- 
termine the junction of the current segment with an adjacent seg- 
ment. In a fast trajectory, for example, the calibrated airspeed V cas is 
monitored during the backward integration of the level-deceleration 
segment. When V cas reaches a certain descent V cas value, integration 
of the level-deceleration segment ends and that of the constant V C3S 
descent segment begins. All other segments are connected similarly. 
Therefore, a trajectory calculation is successful only if the monitored 
variables can reach or capture their desired values. 

In a given descent profile, each segment normally has one vari- 
able that needs to be monitored. These variables are called capture 
variables. When the same segment is put into a different profile, it 
may require a different capture variable. Table 1 lists all possible 
capture variables for each segment. A capture variable is chosen 
from Mach, calibrated airspeed, altitude, or path distance. 

Segment connections will be successful if the capture variable in 
a segment changes monotonically. Let x represent the state vector 
of a segment consisting of some or all of { VJ, h, 5}. The capture 
variable is y and is thus one of {M, V cas , h, j}. Then, equations of a 
segment can be described in the following compact form: 

* = /(*i0 ti < t < r* x(tj) known (9) 

y = v(x) (io) 

The time can go either forward or backward. The final time is defined 
through 


Derivation of Capture Conditions 

In this section, we derive all possible capture conditions for every 
segment. Physical meanings of these capture conditions are dis- 
cussed, and alternative flight strategies are suggested when the con- 
ditions are not met. 

Path distance capture is discussed first. Path distance s is calcu- 
lated from the numerical integration of Eq. (3). The square root in 
Eq. (3) must be well defined, and the ground speed must be posi- 
tive to capture a specified path distance. These requirements impose 
different limits on headwind and tailwind. It turns out that if there 
is a positive constant r < 1 such that 

Vw < rV t (14) 

the square root will be meaningful and 

0 < (1 - r)V t <V g <(\ + r)V r (15) 

A wind magnitude smaller than the true airspeed guarantees path 
distance capture. Path distance capture is needed for every segment. 
Therefore, Eq. (14) must be satisfied in every segment. In most 
flight cases, this condition is satisfied. When predicted wind speed 
becomes comparable to the true airspeed, the airplane should not 
attempt to fly through the region, for safety reasons. 

A. Cruise Segment 

For the cruise segment, Eq. (3) is integrated either forward or 
backward to capture a specified path distance value. Altitude and 
speed (Mach or V^) are held constant. Therefore, the capture condi- 
tion in a cruise segment is the same as that stated above for capturing 
the path distance. 

B. Descending-Acceleration Segment 

In a descending-acceleration segment, true airspeed, altitude, and 
path distance are all integrated numerically. Engine thrust is set at 
minimum (idle) and wind gradient is ignored. True airspeed and 
path distance are integrated from Eqs. (1) and (3), respectively, and 
altitude history is integrated by either specifying a desired descent 
rate h or from Eq. (2) with a desired inertial flight-path angle y*. 
Mathematically, 

h=h* or y*V K (16) 

Correspondingly, 

A* y*V s 

= or ir 07) 

Possible capture variables are Mach number, calibrated airspeed, 
and altitude (Table 1). 

Capturing a Mach number requires M > S m . Differentiating 
Eq. (6) gives 


y(t*) = y* (11) 

Numerically, a capture is declared if 

\y(t*)-y*\<€ (12) 

For the ease of discussion, we assume y(t { ) < y*. Then, a sufficient 
capture condition requires 


m = — tz £ (18) 

VlAR&W) 2 VTa R&( h)\ © ) U ’ 

By combining Eqs. (1), (2), and (18), we have 

dO/dAN] 

0 )\ 

+ m& m yf\AR®(h) ( 19 ) 


T(h, M\ Ka ) > D(h , V,-, L) + mh 


\ L + ~( 
L V, 2 ( 


y > S > 0 (13) 

Monotonicities of capture variables are sufficient but not neces- 
sary for segment connections. A capture variable could vary in an 
arbitrary manner and still reach a specified value before an allowed 
time. However, it is highly desirable to have predictable capture 
variables both for trajectory-synthesis algorithms and in practical 
flights. 


Capturing a calibrated airspeed requires V cxs > S e . From Eqs. (6) 
and (7), we have 
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" 1>0 <21) 

foAh. V l )±(&/p)(V at /V,)(yt/lR& + 1)’ 25 

( pz + 1)^ > 0 (22) 

Therefore 

T(h, M\ /c„) > D(h, V,; L) + mS c f DA (h, V,) 

+m *[‘K + y(^^) +£(Vf ' ft) . (23) 

In a descent flight, h < 0. For subsonic airspeed and standard 

atmosphere, 


Ya = Y a . 


v; 01 


j L+ v l (mw\ >0 
V, 2 V © ) 


E(V„h)± 


( 1 +z)^ 


m RG(h)~ 

dp/dh 

V, 

P(h) _ 


Therefore, conditions in Eqs. (19) and (23) require that the minimum 
thrust in a descending-acceleration flight be larger than a quantity 
that is smaller than the aerodynamic drag. Physically, there must 
be sufficient excess energy for an aircraft to increase Mach number 
or calibrated airspeed. When the minimum thrust is used during the 
descending flight, the aircraft extracts potential energy to gain Mach 
number or calibrated speed. 

If Eq. (19) is not satisfied during the trajectory synthesis, mini- 
mum thrust with specified descent rate h does not provide sufficient 
energy for the aircraft to increase Mach number. One can specify 
a larger rate of descent or a higher thrust setting. It is also possible 
to specify a desired calibrated airspeed for capture. The minimum 
thrust required in Eq. (19) is larger than that in Eq. (23) . As a result, 
it is more difficult to capture a Mach number than a calibrated air- 
speed during a descent acceleration. If Eq. (23) is not satisfied, one 
should either specify a larger rate of descent or use a higher thrust 
setting. Under the same conditions, a heavier aircraft is more likely 
to satisfy Eq.(19) or (23) than a lighter aircraft. 

Capturing a specified altitude is automatic if a desired descent 
rate is selected so that h* < -8 h . If a desired inertial flight-path 
angle is selected, then y* V g < -<5*, which results in 

V s > -- 7 - (26) 

\Y-\ 

This is the same condition as in path distance capture in Eq. (15), if 
8 h is properly selected. It requires that wind magnitude be smaller 
than the true airspeed, as in Eq. (14). In capturing a specified alti- 
tude, Eq. (19) or (23) still needs to be satisfied for Mach number or 
V cas to increase. The desired value of descent rate or inertial flight- 
path angle affects the specific requirement on thrust minus drag in 
Eq. (19) or (23). 

C. Constant Mach Descent Segment 

In a constant Mach descent segment, altitude and path distance 
are numerically integrated and true airspeed is solved as a function 
of altitude from Eq. ( 6 ). We have 

(27) 

Altitude is integrated with one of the following three choices: 1) 
setting the engine thrust at minimum, 2 ) following a desired de- 
scent rate, and 3 ) following a desired inertial flight-path angle. 
Mathematically, 


where y*(h\ Mo) corresponds to setting engine thrust at minimum 
and is given by 


... ... (j\\ mMo-.g n )-p(ft 

Ya (A. Mn) I I { VKdv/dh) _ (d/dA)( v^ 


t. 1 ) ) 

COS 0rw)] + g J 


Possible capture variables are calibrated airspeed and altitude. 

Differentiating Eq. (7) while keeping Mach number constant 
gives 


h /dp\ 
! f MD {h,M (l )\dh) 


f MD (h,MA±-^=[pWB+ 1 ]* 

u v A Wj 

X yjlp(h)B + 1]3 - 1 > 0 (32) 

B = (l + 0.2 Mq ) 3 5 — 1 > 0 (33) 

Capturing a desired calibrated airspeed requires > 8 cy which 
gives 

dp < < o 

d h - c h 

where the descent rate h is given by Eq. (28). In other words, mea- 
sured pressures must decrease as altitude increases. 

Equation (34) defines a bound for the altitude gradient of mea- 
sured pressure at each altitude level. The specific bound depends 
on which one of h y or y, is set to constant. Because atmospheric 
pressure should decrease with altitude, this equation is normally 
satisfied. Actually, it can be used to smooth measurement data of 
atmospheric pressure for trajectory integrations. 

Altitude capture requires h < — The three choices for integra- 
tion of altitude are discussed next. 

1) If the engine thrust setting is held constant, we have 

h = y*V t < -8 h (35) 

which requires that y*(ft; Mo) be negative. In other words, the nu- 
merator and the denominator in Eq. (30) should have different signs. 
In most flight conditions, the minimum thrust is less than the drag: 


T{h % M\ fCo) — D(h y V,; L) <0 
Therefore it is required that 

+ s >0 

L dh dh 

From Eqs. (30), (35), and (36), we obtain 


A ( v„ CO s(W<£ + ^- = S.(A) (37) 

d/i Vt cm 

d g dV t 

— (V„ cos #„,)>- + — 

+ r(/i - ~ £>(fe ' V,; L) = Bj(h) (38) 

mi h 

In addition, a mathematical singularity occurs in the calculation of 
altitude rate if 

■jr(V u , cos 6 ™) -»■ B| (h) (39) 

dh 
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Fig. 2 Wind gradient effect on thrust minus drag. 


In a constant descent segment, altitude and path distance are 
numerically integrated and true airspeed is solved as a function of 
altitude from Eqs. (6) and (8) with as a constant. Define 






- 1 > 0 


We have 


V,(h) = yJlRQ(h)yJ{A/p(h) 


+1)7-1 


(42) 


(43) 


Airspeed and its gradient are functions of altitude for a given Mach 
number; thus the right-hand sides of Eqs. (37) and (38) are functions 
of altitude and atmospheric temperature. The last term on the right- 
hand side of Eq. (38) is negative. These equations establish a feasible 
wind gradient range as a function of altitude for known temperature 
profiles. These limits are satisfied by small wind magnitudes. The 
case of a minimum thrust larger than the drag can be similarly 
studied. 

If these limits are violated, altitude capture in the trajectory syn- 
thesis may fail. Physically, the aircraft has to climb or descend slowly 
to maintain a constant Mach number by taking advantage of the 
wind gradient. If the singularity in Eq. (39) occurs, the calculated 
altitude rate will be very large. There are two alternative methods 
to proceed with trajectory calculations, depending on the nature of 
wind-gradient-limit violations. 

If measured wind gradients violate Eq. (37) or (38) only for a short 
period of time, wind data can be smoothed to satisfy these equations. 
If the violations occur over a long period of time, the aircraft may 
lose its physical ability to follow a constant Mach number while 
descending with minimum thrust. In either case, a desired descent 
rate or inertial flight-path angle can be followed. Correspondingly, 
the required thrust is higher than the minimum if Eq. (37) is vio- 
lated, and a larger drag is needed if Eq. (38) is violated. Details are 
discussed below. 

2) When a descent rate A * is to be followed, altitude capture is 
guaranteed if 


and 


IS 

1 

( d©/d h \ 

| RQ(h)A | 

( dp/dh \ 

d h 2 ' 

l © J 

V',[A/p(A) + l]f * 

V P 2 ) 


Altitude is integrated from Eq. (28). The corresponding air- 
relative flight-path angle is given in Eq. (29). Possible capture vari- 
ables are Mach number and altitude. In this case, y* is given by 

* {k . y x 1 [ T(h,M-,K 0 )- D(h, V,;L) 

° ’ C “ m\ V t [(dV'/dh)-(d/dh)(V w cose m )] + g 


(45) 

where the Mach number changes according to Eq. (8). Actually, we 
have, from Eqs. (6) and (43), 


M = ^ y j\AJp(h)J r l]f - 1 


(46) 


and 


A 

' dp/dh " 

_M(A/p + 1)$_ 

P 2 


(47) 


A* < 


(40) 


Capturing a desired Mach number requires M < which 
results in 


3) If a desired inertial flight-path angle is to be followed, altitude 
capture requires 

h = y-V g < -S h =*■ V s > (S h / - y’) (41) 

which is the same condition as path distance capture in Eq. (15), 
if the 5* is properly selected. The wind magnitude must be smaller 
than the true airspeed. 

Figure 2 shows required ranges of (T—D) in various cases of wind 
gradients to ensure a descending flight with a constant Mach number. 
If Eq. (37) is violated, a descending aircraft would experience a large 
decrease in headwind. The aircraft would lose airspeed if it keeps 
descending with minimum thrust. Mathematically, the aircraft could 
climb to take advantage of the increasing headwind. In an actual 
flight, the pilot has to use a larger thrust than idle to maintain a 
constant Mach number. In contrast, a larger drag is needed if Eq. 
(38) is violated. For example, spoilers can be used to augment the 
aerodynamic drag. 

Following a specified descent rate or inertial flight-path angle 
requires a thrust level different from the minimum thrust. If both 
Eqs. (37) and (38) are satisfied, the computed thrust is less than the 
drag. When either Eq. (37) or (38) is violated, the computed thrust 
minus drag follows the trends in Fig. 2. In particular, the required 
thrust is close to drag if the measured wind gradient is close to the 
upper bound as in Eq. (39). 

D. Constant Calibrated Airspeed Descent Segment 

A constant V cas descent segment is similar to a constant Mach 
descent segment, except that the calibrated airspeed is maintained 
constant instead of the Mach number. As a result, numerical inte- 
grations capture a desired Mach instead of a desired V cas . 


or 


dp „ . p 2 fco(h ; V CM ) „ 

dh~ Sm h <0 




(48) 


(49) 


where h is given in Eq. (28) and 

fco(h ; V caJ ) = 0/A*JS)(A/p+\)$,J (A/p + 1)^-1 > 0 (50) 


Again, the altitude gradient of measured pressure profile must be 
negative and stay within a certain bound. The specific bound de- 
pends on which one of a. A, or y t is set to constant. If this bound 
is violated, smoothing of the measured pressure data to continue 
trajectory synthesis is recommended. 

Altitude capture requires A < -<5* , which results in the same con- 
ditions as in the constant Mach segment, that is, Eqs. (37) and (38) 
for constant engine thrust setting, Eq. (40) for following a constant 
descent rate, and Eq. (41) for following a desired inertial flight-path 
angle. Alternative strategies and thrust-minus -drag requirements are 
similar as in the case of a constant Mach descent. However, the true 
airspeed gradient in Eq. (44) for a constant descent is differ- 
ent from that in Eq. (27) for a constant Mach descent. As a result, 
Eqs. (37) and (38) give different bounds in the two cases. Although 
details are omitted, capture requirements for a constant descent 
allow for larger wind gradients than those in a constant Mach descent 
at the same airspeed and altitude. 
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E. Level-Deceleration Segment 

In a level -deceleration segment, true airspeed and path distance 
are integrated from Eqs. ( 1 ) and (3) and altitude is held constant. En- 
gine thrust is set at minimum. A level -deceleration segment can be 
connected to another segment through the capture of either path dis- 
tance or speed (Table I). Usually, Mach number is used for capture 
at high altitudes whereas calibrated airspeed is used at low altitudes. 

Capturing a calibrated airspeed in a level deceleration segment 
requires V cas < — S c . From Eqs. (6) and (7), we have 


Kas = 




/ ld (^ 0, Vr) 


(51) 


where 



As a result, 

T(h 0 , M ; K {) ) < D(h 0 , V t ; L) - mS c f LD (k 0 , V,) (53) 

Capturing a Mach number in a level -dec deration segment re- 
quires M < —& m - From Eq. (6), we obtain 

T{h 0y M\Ko)< D(h 0 , V f ; L) - mS m ^\ARQ 0 (54) 

Therefore the minimum thrust must be smaller than the aerody- 
namic drag by a certain margin in a level-deceleration flight. This 
margin is different in capturing a Mach number from that in cap- 
turing a calibrated airspeed. This requirement is not in conflict with 
Eqs. (19) and (23) for the descent-acceleration segment. Actually, 
all these requirements together define a feasible interval for the min- 
imum engine thrust. If Eq. (53) or (54) is not satisfied, a larger drag 
needs to be used. 


conditions apply directly if the integration time steps are reasonably 
small. Furthermore, the continuous-time conditions with larger val- 
ues of S will guarantee captures in numerical integrations. 

The discrete capture conditions are more complicated than the 
continuous-time conditions, making them difficult to use. Therefore, 
the continuous-time conditions are recommended for use in these 
applications. 

Applications of Capture Conditions 

The capture conditions have direct applications in trajectory syn- 
theses for air traffic control automation. They can be used in the 
trajectory synthesis process to verify convergence conditions and to 
suggest alternatives in case a capture fails. They also can be used 
to smooth measured data of atmospheric conditions. In these two 
applications, values of S need to be specified. Specific values of S de- 
pend on time and iteration allowed for convergence. The following 
expression gives a useful guide: 


Nmax A^min 

where Nmax is the maximum number of integration steps allowed, 
and A t mitt is the smallest practical integration step size. 

In addition, capture conditions can be used for analyzing feasible 
ranges of flight altitudes and speeds both for actual flight and for 
trajectory synthesis. In this case, S can be set to zero for conve- 
nience. Equations (1-3) represent approximate point-mass aircraft 
equations. Derivations of capture conditions are similar even if more 
accurate equations are used. 

The specific limits imposed by the above capture conditions de- 
pend on the thrust and aerodynamic models. Therefore, these limits 
vary from one aircraft to another. On the other hand, the derived 
capture conditions represent physical requirements on aircraft per- 
formance and atmospheric environment. These conditions are valid, 
independent of specific algorithms used for trajectory synthesis. 

The above-derived capture conditions also can be used if two cap- 
ture variables are simultaneously searched in segment connections. 
Capture conditions for both variables must be satisfied. 


Effects of Numerical Integration Scheme 

In actual implementations of trajectory synthesis, numerical in- 
tegrations are performed. The current CTAS uses a second-order 
Runge— Kutta scheme. 8 Based on Eqs. (9) and(10), the numerical 
integrations calculate a sequence of states from 

Xj + i = Xj + AtF(Xj, tj\ At) (55) 

y; + ! = y(x j + i) (56) 

for j = 1 N, where 

F{xj , ty. At) = ±{/(x,, t } ) + f[xj + Atf(xj, tj ), f, + 1 ]} (57) 

The capture logic in CTAS at each integration step follows. Again, 
we assume y(fi) < y* for the convenience of discussion. 

1) Start with a At and use Eqs. (55) through (57) to calculate 
values of state and capture variables at the next step. 

2) If \y(t j + i) - y*l < e, the variable is captured and t* = t j+ j. 
The integration of the current segment is completed. 

3) Otherwise, if 

y*-y0j) >€ and y(t j + \)-y*>€ 

then. At is too large. Decrease A t, reassign initial conditions, and 
go to step 1. 

4) Or, assign tj + At tj and advance one step. 

The current capture accuracies in CTAS are 0.0003 for Mach 
number, 0.1 kn for calibrated airspeed, 20 ft for altitude, and 0.005 
n mile for path distance. 

Capture requirements derived for Eqs. (55-57) are different 
from conditions derived above using continuous-time derivatives. 
Nonetheless, a segment-by-segment analysis shows that the discrete 
conditions converge to the continuous-time conditions as the inte- 
gration time step shrinks to zero. In other words, the continuous-time 


Profile Capture Conditions and Example 

Conditions are derived for every segment for all possible capture 
variables. The defined descent profiles, such as the fast type illus- 
trated in Fig. 1, combine a series of segments in a predetermined 
order. The capture variable for each segment in these profiles is usu- 
ally fixed. Therefore, one can group capture conditions for all the 
segments in a profile to obtain profile capture conditions. 

For the fast-type descent profile in Fig. 1 , the level -deceleration 
segment captures a descent in backward integration, the constant 

V cas descent segment captures the descent Mach number in backward 
integration, the constant Mach descent segment captures the initial 
altitude in backward integration, and the descending -acceleration 
segment captures the descent Mach number in forward integration. 
Finally, the trajectory calculation is completed by a backward in- 
tegration of the cruise segment to capture the initial path distance. 
Other profiles can be analyzed similarly. 

We now summarize the capture conditions for a jet airplane 
flying a fast-descent-profile trajectory. The thrust for a jet engine 
is controlled by using the engine pressure ratio, which is set at 
minimum (idle) during the descending-acceleration segment, the 
constant Mach descent, the constant V cas descent, and the level- 
deceleration segment. During the descending-acceleration flight, a 
desired descent rate is followed. Capture requirements with 5 = 0 
are as follows: 

1) The wind magnitude should always be smaller than the true 
airspeed as required in Eq. (14). 

2) During the descending-acceleration segment, the minimum 
thrust needs to satisfy Eq. ( 1 9) to capture a specified Mach number or 
Eq. (23) to capture a specified calibrated airspeed. Mathematically, 
capturing a Mach requires 


T(h, M, *o) - D(h, V,-, L) J g ,V'i 

( d0/d/t V 

m [v, + 2 ' 

l © )\ 
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Fig. 3 Capture conditions in descending-acceleration segment 


and capturing a calibrated airspeed requires 
T (h, M\ xrp) - D(h, V t ; L) 
m 





>0 


(60) 


where E(V n h) > 0 is defined in Eq. (25). 

3) During the constant Mach descent, the wind gradient must stay 
below the limit defined by Eq. (37) to capture a specified altitude. 
Otherwise, a higher thrust is required. 

4) During the constant V cas descent, the wind gradient must stay 
below the limit defined by Eq. (37) to capture a specified altitude. 
Otherwise, a higher thrust is required. To capture a specified Mach 
number, the altitude gradient of measured pressure needs to be neg- 
ative, as required by Eq. (48) or (49). 

5) During the level-deceleration segment, the minimum thrust 
needs to be less than the drag as in Eq. (53) to capture a specified 
calibrated airspeed. 

The above conditions are now illustrated using a Boeing 727-200 
aircraft in descent. Powered by three Pratt and Whitney JT8D-15 
engines, this aircraft has a typical descent weight of 140,000 lb and 
a default descent calibrated airspeed of 280 kn. The maximum V cas 
is 350 kn and the minimum V cas is 230 kn. Maximum-cruise Mach 
number is 0.82. In Figs. 3-6, the standard atmospheric pressure, 
density, and temperature are assumed. A descent rate of 3000 ft/min 
is used in segments with specified descent rate. 

In the CTAS aircraft models, 9 drag coefficients for the Boeing 
727-200 aircraft are given in a two-dimensional table at 6 Mach 
numbers and, for each Mach number, 11 lift coefficients. The min- 
imum thrust is given in a two-dimensional table at 10 altitudes and 
10 Mach numbers. In this paper, the drag coefficient is fitted into 
a fifth-order polynomial of lift coefficient at each Mach number. 
A linear interpolation is used for drag coefficient between any two 
specified Mach numbers. A two-dimensional linear interpolation is 
used to determine the minimum thrust between specified values of 
altitude and/or Mach number. 

Figure 3 plots the capture requirements in Eqs. (59) and (60) 
for a descending-acceleration segment at various Mach numbers 
and three different altitudes. Positive regions of the curves in Fig. 3 
define a feasible range of Mach numbers and altitudes within which 
the descending-acceleration integration can capture. For example, 
for h > 40,000 ft, the maximum Mach number that can be captured 
is around 0.814. Clearly, it is easier to capture a calibrated airspeed 
than a Mach number at the same altitude and speed. These curves 
will shift up as the aircraft weight or the descent rate increases. If 
capturing a Mach number becomes difficult, one can either increase 
the descent rate or switch to capture a calibrated airspeed. 

Figures 4 and 5 show the upper bound on wind gradient in Eq. (37) 
for a constant Mach descent and a constant V CiS descent, respectively. 
These figures show different trends for the constant Mach descent 
and for the constant V C25 descent as altitude increases. The maxi- 
mum wind gradient that can be tolerated decreases as flight speed 
increases. Over the altitude range of 10,000-35,000 ft, a wind gra- 
dient less than 0.03 s -1 meets the requirement at all speeds. This is 
consistent with the findings of numerical experiments. 



Fig. 4 Upper bound on wind gradient in constant Mach segment. 



Fig. 5 Upper bound on wind gradient in constant V c „ segment 



Finally, Fig. 6 shows the minimum thrust minus drag at different 
altitudes and speeds. This quantity is always negative in this exam- 
ple. It is always possible to capture a calibrated airspeed within the 
range of 230-350 kn. 

In summary, this example shows that the capture conditions are 
satisfied for a Boeing 727-200 in the standard atmospheric con- 
ditions, except possibly in capturing a Mach number during the 
descending-acceleration segment. 

Conclusions 

This paper presents capture conditions that affect the successful 
connections of flight segments in numerical integrations for tra- 
jectory synthesis. In a commercial descent, an aircraft trajectory 
is divided into a series of segments according to preferred pilot 
procedure. A capture variable, either path distance, Mach number, 
calibrated airspeed, or altitude, is monitored in one segment, and the 
next segment begins when the capture variable reaches a specified 
value. Successful connections of segments in trajectory integrations 
are guaranteed if capture variables are monotonic in the segment for 
which they are monitored. In this paper, monotonicity conditions for 
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capture variables in different segments of various descent profiles 
are derived. Alternative flight strategies are suggested when some 
capture conditions are not satisfied. 

These conditions impose bounds on measurements of atmo- 
spheric wind, pressure, and temperature, as well as thrust minus 
drag. These bounds can be used to smooth measurement data for 
trajectory synthesis. They can be used in trajectory-synthesis algo- 
rithms to provide convergence checks and alternative instructions. 
For given aircraft performance models, these capture conditions 
define feasible ranges of altitudes and airspeeds within Which an 
aircraft can be directed to fly. 

The example of a Boeing 727-200 aircraft descending in a stan- 
dard atmosphere shows that the capture conditions are satisfied in 
most flight segments, except possibly the descending-acceleration 
segment when capturing a Mach number. 
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